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Abstract: 

Relativistic mean field (RMF) theory of nuclear matter with the isovector scalar mean 
field corresponding to the 5-meson [ao(980)] is studied. While the 5-meson mean field 
vanishes in symmetric nuclear matter, it can influence properties of asymmetric nuclear 
matter in neutron stars. The RMF contribution due to 5-field to the nuclear symmetry 
energy is negative. To fit the empirical value, E s w 30MeV, a stronger p-meson coupling 
is required than in the absence of the 5-field. The energy per particle of neutron matter 
is then larger at high densities than the one with no 5-field included. Also, the proton 
fraction of /3-stable matter increases. Splitting of proton and neutron effective masses due 
to the 5-field can affect transport properties of neutron star matter. 
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RMF models of nuclear matter [1] play an important role in dense matter calcula- 
tions, where relativistic effects increase with density. They represent a relativistic Hartree 
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approximation to the one-boson exchange (OBE) theory of nuclear interactions, with ap- 
propriately readjusted couplings to fit the saturation properties of nuclear matter. In the 
OBE description of nucleon-nucleon scattering, isoscalar mesons, a, rj, and u>, and isovector 
mesons, 7r, p, and 5, are involved. 

In the RMF approximation, mean fields of pseudoscalar mesons rj and 7r vanish in 
normal nuclear matter. The standard RMF model of nuclear matter [1] includes isoscalar 
mesons a and u, and the isovector p-meson. However, no 5-meson contribution is present, 
although for asymmetric nuclear matter the density to which the 5-field can couple does 
not vanish, < >^ 0. 

The 5-field is not expected to play an important role for nuclei, whose isospin asym- 
metry is not large. Also, the short range of the 5-meson exchange justifies neglecting its 
contribution at saturation density. However, for strongly isospin-asymmetric matter at 
high densities in neutron stars the contribution of the 5-field should be considered. It 
is our aim to include the 5-meson into the RMF model and study consequences of this 
generalization of the RMF theory. To our knowledge, such an extension of RMF models 
has not been studied. 

We assume here that all the meson fields have the Yukawa couplings to nucleons. The 
interaction lagrangian reads 



Here 6 is the isovector scalar field of the 5-meson. The free field lagrangians for ip, a, u, 
and p fields are the same as in Ref.[l]. For the 5-field we use the standard lagrangian 




(1) 



U 



= \djd»5- 1 -m]$ 2 . 



(2) 



For the cx-field we adopt the potential energy term of Boguta and Bodmer [2], 
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U(a) = ^bma 3 + l -ca\ (3) 

where m is the bare nucleon mass. 

The relevant components of meson fields are cr, ujq, and for the 5-field the isospin 
component 5^ 3 \ All remaining components vanish, in particular 5^ = 5^ = 0. The 
values of mean fields in the ground state are determined by proton and neutron densities. 
The field equations for vector meson fields give u>q = (g^/m^nB and 
(ffi = (9p/ m f>)(^ x ~ l) n B, where jib is the baryon density and x = np/np is the proton 
fraction. The field equations for scalar fields are non-trivial: 



dll 

mla+ — =g a (n s P + n s N ), (4) 

and 



m 2 s m =gs(n s P -n s N ). (5) 
Here np and n s N is, respectively, proton and neutron scalar density: 



{ 2n > Jo ^k 2 + mf 



In Eq.(6) mp and tun is, respectively, proton and neutron effective mass, 



mp=m - 9a a - g s 5 (3 \ 



(7) 



and 



m N = m- g a o + g s S^ . 

The energy density of uniform nucleon matter is 
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(8) 
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The model parameters in the isoscalar sector, C 2 = g^/m^C 2 = <7 2 /m 2 ,& = b/g 3 , 
and c = c/g 3 , are adjusted to fit the saturation properties of symmetric nuclear matter, the 
saturation density no = 0.145/m -3 , the binding energy wo = —16MeV per nucleon, and 
the compressibility modulus Ky 280MeV. The fourth parameter, e.g. c, can be used to 
measure stiffness of the equation of state of symmetric nuclear matter. In the following we 
use two sets of parameters which reproduce the saturation properties but differ at higher 
densities. The soft equation of state is specified by the parameters C 2 = 1.582/m 2 ,C 2 = 
1.019/m 2 ,6 = —0.7188, and c = 6.563. For the stiff equation of state the parameters are 
Cl = 11.25/m 2 , C 2 = 6.483/m 2 , b = 0.003825, and c = 3.5 x 10" 6 . 

In the spirit of RMF models, the parameters C 2 = g 2 p /m 2 p and C 2 = g 2 /m 2 s of the 
isovector sector should be constrained to fit the nuclear symmetry energy, E s = 3l±4MeV 
[3]. In terms of the model parameters, the symmetry energy is 

ft - -r 2 n i k o r 2 m o n ° 

S ^ P ° Q^/kJ^ 2 5 2(k 2 + ml)(l + C]A(k ,mo)y ^ 

where 

A(k ,m ) = ^3 ^ {p2 + mlr/2 (ID 

is a function of the Fermi momentum, ko = kp = /cat, and the effective mass, mo = mp = 
raw, of symmetric nuclear matter at saturation density. The constraint (10) gives C 2 as 
a function of Cf , C 2 = C 2 (Cf ), which is shown in Fig.l. Hence by demanding to fit the 
bulk properties of nuclear matter we are not able to fix independently C 2 and Cf. Instead, 
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we explore here a range of values of Cf corresponding to the Bonn potentials A, B and C 
[4]. The maximum value Cf = 2.6/m 2 corresponds to the parameter gg = 8.0 of the Bonn 
potential C [4]. 

The fact that the constraint (10) gives C 2 as a monotonically increasing function of 
Cf plays a crucial role in our analysis. The nuclear symmetry energy, Eq.(10), is derived 
in the RMF approximation which neglects exchange contributions. Inclusion of the Fock 
terms in relativistic nuclear matter calculations gives a sizable correction to the nuclear 
symmetry energy [5]. We expect, however, that adding the Fock contribution will not 
change the general behaviour of the curves in Fig.l, although it can affect their intercept 
and slope. This topic will be discussed in detail elsewhere [6] . 

To estimate the effect of the 5-field we first consider the case when the coupling con- 
stant C 2 is adjusted with no 5-field present (Cf = 0). For the soft and stiff equation of 
state we find, respectively, C 2 = 5.0/m 2 and C 2 = 4.29/m 2 . In Fig. 2 we show correspond- 
ing energies of pure neutron matter. When the 5-field is switched on, the energy/particle 
becomes smaller. In Fig. 2 we show curves for a few values of Cf from the range compatible 
with the Bonn potentials A, B and C [4]. 

The contribution due to the 5-field is attractive and quite large. In case of the stiff 
equation of state the neutron matter becomes selfbound for Cf > 1.1/m 2 . Already a 
smaller value, Cf « 1.0/m 2 , lowers the symmetry energy to about 20MeV. One should 
note that for C 2 = 2.5 fm 2 the neutron matter becomes as strongly bound as the symmetric 
nuclear matter. For the soft equation of state the binding is stronger. 

To avoid such an unphysical situation, the increased binding due to the 5-field has 
to be balanced by the higher repulsion due to the p-field. In particular, if we demand 
the symmetry energy to be reproduced, the formula (10) shows that the parameter C 2 
has its lowest value for Cf = 0. For any finite value of the 5-coupling, Cf > 0, the 
strength of the p-coupling, C 2 , increases. In this case inclusion of the 5-meson results in 
higher energy /particle at high densities, where the contribution of vector mesons dom- 
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inates. The value C 2 (0) should be regarded as the lower bound, with the actual value 
larger, corresponding to some finite Cf . Similarly, the energy/particle of neutron matter 
is bounded from below by values obtained for C 2 (0). Actual energies, for finite Cf, are 
higher. Below we show results for Cf = 2.5/m 2 which is close to the value of the ratio 
gj/m'g corresponding to the Bonn potential C [4]. 

The presence of the mean 5-field leads to splitting of proton and neutron effective 
masses, Eqs.(7) and (8). In Fig. 3 and Fig. 4 we show the effective masses mp and m^ 
as functions of baryon density for a few values of the proton fraction x. It is interesting 
to note that for the stiff equation of state, Fig. 3, the proton and neutron effective masses 
in pure neutron matter both decrease monotonically with density. For the soft equation 
of state, Fig. 4, the proton effective mass in pure neutron matter increases with density, 
approaching asymptotically the value mp = 2m for the largest couplings Cf . In contrast, 
the neutron effective mass decreases monotonically. Such a splitting of proton and neutron 
mass can affect the transport properties of dense matter. 

The energy per particle is presented in Fig. 5 where we show results for pure neutron 
matter for both soft and stiff equations of state. For comparison, curves for neutron matter 
with no 5-field included, and for symmetric nuclear nuclear matter, are also shown. One 
can notice that the energy /particle with the 5-contribution included increases more rapidly 
with neutron matter density for both equations of state than the energy for Cf = 0. 

The proton fraction of ^-stable neutron star matter, which satisfies the condition 

m ~ = yj 'k 2 N + m% - yj 'hp + mp + ^C 2 (l - 2x)n B = A*e, (12) 

where the electron chemical potential is /j e = (S^xup) 1 ^ 3 , is shown in Fig. 6. One can 
notice that for Cf = 2.5/m 2 the proton fraction at high densities is larger than for Cf = 
for both equations of state. It exceeds the critical value for the direct URCA process to 
be of importance in the cooling of neutron stars, which is xjjrca ~ 0.11. 

In conclusion, the RMF theory with the contribution due to the 5-meson mean field, 
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constrained to fit bulk properties of nuclear matter, predicts higher energy /particle of 
the neutron matter at a given baryon density than in absence of the 5-field. Although 
the 5-field provides an additional binding, fitting of the symmetry energy, E s « 30MeV, 
requires a stronger p-meson coupling which dominates at higher densities. The proton 
fraction of the ^-stable neutron star matter is considerably higher than the one with no 
5-field contribution. Proton and neutron effective masses are split, an effect which can 
modify the transport properties of dense matter. 

This research is partially supported by the Polish State Committee for Scientific Re- 
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References 

[1] B. D. Serot and J. D. Walecka, Adv. Nucl. Phys. 16 (1986) 1. 

[2] J. Boguta and A. Bodmer, Nucl. Phys. A292 (1977) 413. 

[3] W. D. Myers and W. D. Swiatecki, Ann. Phys. 84 (1973) 186; 

J. M. Pearson, Y. Aboussir, A. K. Dutta, R. C. Nayak, M. Farine, and F. Tondeur, 

Nucl. Phys. A528 (1991) 1; 

P. Moller and J. R. Nix, At. Data and Nucl. Data Tables 39 (1988) 219; 

P. Moller, W. D. Myers, W. J. Swiatecki, and J. Treiner, At. Data and Nucl. Data 

Tables 39 (1988) 225; 

W. D. Myers, W. J. Swiatecki, T. Kodama, L. J. El-Jaick, and E. R. Hilf, 
Phys. Rev. C15 (1977) 2032. 
[4] R. Machleidt, Adv. Nucl. Phys. 19 (1989) 189. 

[5] A. Bouyssy, J.-F. Mathiot, Nguyen Van Giai, and S. Marcos, Phys. Rev. C 36, 
(1987) 380; 

M. Lopez-Quelle, S. Marcos, R. Niembro, A. Bouyssy, and Nguyen Van Giai, Nucl. 
Phys. A483 (1988) 479. 
[6] S. Kubis and M. Kutschera, in preparation. 

7 




8 



Fig.2 



200 



150 - 



> 

CD 



400 - 



TO 



50 - 








■50 







c 2 



C 2 = 2.5 fm 2 



_i_ 



2 

nB/n 



The energy per particle of pure neutron matter for a few values of the 5-meson coupling, 
Cf . For the same value is used for all curves; it is the value which fits the nuclear 
symmetry energy with no 5-field. Solid and dashed curves correspond to the soft and 
stiff equation of state, respectively. For comparison, the energy per particle of symmetric 
nuclear matter is shown (dash-dotted curves). 
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Effective masses of protons and neutrons for a few values of the proton fraction, for the 
stiff equation of state and for C| = 2.5/m 2 . 
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Fig.4 



The same as Fig. 3 for the soft equation of state. 
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Fig. 5 
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The energy per particle of neutron matter for Cf = 2.5/m 2 , corresponding to the Bonn 
potential C, and with no 5-meson contribution. The parameter C 2 p is adjusted to fit the 
nuclear symmetry energy in both cases. Solid and dotted curves are, respectively, for the 
soft and stiff equation of state. The energy per particle of symmetric nuclear matter is 
also shown. 
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Proton fraction of neutron star matter in the presence of the 5-field and with no 5-field. 
Solid and dotted curves correspond to the soft and stiff equation of state, respectively. 
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